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We describe several new ways of specifying the behaviour of Lemaitre-Tolman (LT) models, in 
each case presenting the method for obtaining the LT arbitrary functions from the given data, and 
the conditions for existence of such solutions. In addition to our previously considered 'boundary 
conditions', the new ones include: a simultaneous big bang, a homogeneous density or velocity 
distribution in the asymptotic future, a simultaneous big crunch, a simultaneous time of maximal 
expansion, a chosen density or velocity distribution in the asymptotic future, only growing or only 
decaying fluctuations. Since these conditions are combined in pairs to specify a particular model, 
this considerably increases the possible ways of designing LT models with desired properties. 
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I. AIM AND MOTIVATION 

The Lemaitre-Tolman (LT) metric is the most widely 
used model of cosmic inhomogeneity, being suitable for 
both large and small scale inhomogeneities. As the sim- 
plest inhomogeneous solution of Einstein's equations, it 
is relatively easy to work with, but retains the full non- 
linearity of the field equations. 

However, because of this nonlinearity, generating mod- 
els with a specific desired evolution was not all that easy, 
and involved a large measure of guesswork in chosing the 
arbitrary functions of the model, as well as repeated nu- 
merical evolution to check if satisfactory results had been 
obtained. 

The usefulness of the LT metric is greatly increased if 
specific models can be designed to have certain prop- 
erties that satisfy observational or theoretical require- 
ments. Previous papers 0,1113, (hereafter papers I, 
II and III) have shown how to generate LT models that 
evolve from a given (spherically symmetric) density or 
velocity profile at time ti to a second given density or 
velocity profile at time i2- 

We here extend the methods by which models may 
be constructed. We seek to solve for the arbitrary func- 
tions E{M) and tsiM) that determine the LT model 
that evolves according to the following requirements, in 
various combinations, where the first two were considered 
previously: 

• a density profile Pi{M) is given at time ti, 
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• a velocity profile {R^t)i{M) is given at time ti, 

• the bang time is simultaneous, 

• the cruch time is simultaneous, 

• the time of maximum expansion is simultaneous, 

• the model becomes homogeneous at late times, 

• only growing modes are present, 

• only decaying modes are present, 

• a velocity profile (i?_f)(M) is given at late times, 

• a time-scaled density profile t^p{M) is given at late 
times. 

It requires two of the above conditions to specify an LT 
model, and in various contexts different combinations 
may be useful. We work them out here in a system- 
atic way for future reference. Where two conditions are 
specified on the same time slice, this often requires spe- 
cial treatment. (In contrast, Q shows how to determine 
the LT functions from observational data on the past null 
cone.) 



II. OUTLINE OF THE LEMAITRE-TOLMAN 
MODEL 

The Lemaitre-Tolman (LT) model is a spherically 
symmetric, nonstatic solution of the Einstein equations 
that is inhomogeneous in the radial direction. The matter 
source is a perfect fiuid with zero pressure, i.e. dust, and 
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the coordinates are comoving with the matter particles. 
The metric is 



1 + 2E{r 



(2.1) 



where E{r) is an arbitrary function that determines the 
local curvature of constant t slices, and R_r is the r deriva- 
tive of the areal radius R{t, r). 

With A = assumed, the evolution equation for R is 



R.t' = 



2M 

'r~ 



2E. 



(2.2) 



where M(r) is a second arbitrary function, representing 
the total gravitational mass within the comoving matter 
shell at constant r. Note that in this equation E{r) has 
a second interpretation as the local energy per unit mass 
of the dust particles, and thus it also determines the type 
of evolution (see below). The density is 



2M^r 
R^Rr 



where n ■ 



SttG 



(2.3) 



The evolution equation (|2.2|) has three types of solution: 
Elliptic, E <0: 



R{t,r) 



M 



(1 - cos 77), 



i-m 

rj — sin ?7 
Parabolic, = 0: 

"Q 

R{t,r) 

Hyperbolic, i? > 0: 
M 



3/2 



M 



M{t^tB{r)f 



-{t^tB{r)); (2.4) 



1/3 



(2.5) 



R{t,r) 



2E 



(cosh?7 — 1), 
sinh77-77= ^-J^^ {t-tB{r)); (2.6) 



where 77 is a parameter and tsir) is the third arbitrary 
function, representing the local time at which the big 



bang occurs. The parabolic evolution is the E ^ 
limit of the other two evolutions, obtained by noting that 
r]/\/—2E remains finite. It is perfectly possible to have 
adjacent elliptic and hyperbolic regions in one model, the 
evolution being parabolic on the boundary where E = 0, 
but in general E' ^ 0. 

Elliptic models have both a big bang t = tB and a big 
crunch. 



tc{r)=tB+T{r), 

where the lifetime of each worldline is 

27rM 
- (-2£;)3/2- 

The instant of maximum expansion is 



tMX (r) ^tB + 



ttM 



(„2£;)3/2' 

at which moment the maximum areal radius 

^ , , M 
RMx(r) = 

is reached. 

The homogeneous case is obtained by setting 
E cx Af , tB = constant 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



and it is a Friedmann model, i.e. a Robertson- Walker 
(RW) model with zero pressure. 
In the following, we will use the notation 



a = R/M^^^, 
X = {-2E)/M^/^, 
X = {2E)/M^/^, 
b = R^t/M^/^. 



E <0, 
E>0, 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



The parametric solutions H2.4|) and H2.6|l may be written, 
for the expanding hyperbolic (HX), expanding elliptic 
(EX), collapsing elliptic (EC), and collapsing hyperbolic 
(HC) cases as 



HX: t = tB + x-^/'^ 1^(1 + xa)2 - 1 - arcosh(l + a;a)| , E > 0, t > ts; (2.16) 

EX: t = tB + x^^l'^ |arccos(l + xa) - ^\ - (1 + xa)^ |, £:<0, 0<77<7r; (2.17) 

EC: t = tB + x^^l"^ |27r - arccos(l + xa) + ^1 - (1 + ^af- | , S < 0, tt < 77 < 27r; (2.18) 

HC: t = tc- x"^/^ I v/(l + xay - 1 - arcosh(l + a;a) j , E > 0, t<tc. (2.19) 
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Alternatively, from 12. 2() and p. 13(1 or H2.14|l we have 6^ = 2/a ± x, so substituting for a in the above gives 
HX: t = tB+ x-^l"^ 

EX: t = tB+x' 



EC: t^tB+x < 27r - arccos 



HC: t = tc-x' 




1 — arcosh 



(2.20) 
(2.21) 
(2.22) 
(2.23) 



r 



There are also two borderline cases, the parabolic case, 
and the elliptic case which has reached maximum ex- 
pansion at time t. Since the above expressions are not 
numerically well-behaved near these borderlines, we will 
present series solutions for near-parabolic (nP), and for 
near-maximum-expansion (nMX) models. 

The time reversed solutions, in which the hyperbolic 
and parabolic solutions are collapsing (HC & nPC), are 
also possible, but not as relevant for cosmology. Thus 
they will only be given for certain cases to make a com- 
plete listing of possibilities. 

It is convenient in what follows to use M(r) as the 
radial coordinate (i.e. r — M{r)), since we are not re- 
ally considering vacuum, or extrema in the spatial sec- 
tion ( "bellies" and "necks" ) . Thus we can integrate (|2.3|) 
along a constant time slice, t ^ ti, to obtain 



Rf =ajM - 



M 
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■AM. 



npi{M) 



(2.24) 



This equation tells us that, if we have a density profile 
Pi{M) > as a function of mass given at a particu- 
lar time ti^ then a straightforward integration gives us 
Ri{M) on that time slice. Ways of coping with regions 
of zero density were discussed in . 

Ideally, we seek LT models that have regular origins, 
and are free of shell crossings and surface layers 0|. 

See Hi as well as 111 la, lal for more details on LT 
models. 



A. Profile to Profile Solutions 

As background to the present work, we very briefly 
summarise the results of the previous papers. 

Paper I showed that, if a spherically symmetric density 
profile is given at two different times, i.e. 

p = pi{M)>Q at t^h, (2.25) 
p = p2(M)>0 at t = t2, (2.26) 



then there always exists a LT model that evolves from 
one to the other. The formulas for the arbitrary functions 
E{M) and tB{M) were given, as well as the conditions 
that determine which type of evolution applies at each 
M value. Since the formula for E{M) can only be given 
implicitly, a numerical code was written to implement 
this algorithm. An example of the formation of an Abell 
cluster with a realistic density profile at the present day, 
starting from a small fluctuation at recombination, was 
calculated and its evolution plotted. 

Although the formulas given were for t2 > ti and for an 
expanding model at ti, it is easy to adapt to the time- 
reverse scenario. A key step in the solution process is 
converting a given density profile Pi{M) to an areal radius 
profile Ri{M) via equation Thus, if Ri{M) were 

given instead, this is easily incorporated into the method. 

Since the two density profiles are arbitrary, it is en- 
tirely possible that the resulting LT model could develop 
shell crossings somewhere, but it is easy to check for this 
once E{M) and tB{M) are known. A second possibility 
to check for is whether 2E has reached —1, which indi- 
cates a maximum in the spatial sections, R_r = has 
been reached. Assuming both given densities are finite 
and non-zero at this point, it would be a regular comov- 
ing maximum, i?.r(i, Afmax) = for all t. Beyond this 
locus, a regular model would have M and R decreas- 
ing with increasing radial distance, so, in order to make 
further progress, one would have to replace p{M) with 
p(2Af,„ax — M). But this eventuality was not included 
in the numerics, as models of objects that large were not 
contemplated. 

Paper II considered the possibility that velocity profiles 
might be given, 



R,t ^ iR,t)^{M) > at t = U 



(2.27) 



and again showed how to find the LT model that evolves 
between two such profiles or between one density and one 
velocity profile. Several numerical examples were given, 
illustrating different possibilities, including an improved 
Abell cluster model and others that showed how radical 



4 



changes in the density profile are possible, and highlight 
the fact that both density and velocity fluctuations at 
recombination play a significant role. There was also a 
model of the development of a void, but the density and 
velocity fluctuations at recombination could not both be 
made small enough in any of the models tried. Improved 
void models were constructed in p^ . using a realistic 
present-day density profile based on observations, and 
various density or velocity profiles at last scattering con- 
sistent with CMB observations, as well as a variety of 
model parameters, particularly rimattcr and fl^, Hq, k. 
It was found that in each model there was always some 
inconsistency with observations: either the void density 
was not low enough, or the initial velocity was too high, 
or the density profile at the void wall was too steep and 
developed a shell crossing. A much improved consistency 
with observations was achieved by Bolejko [T]| by includ- 
ing an inhomogeneous radiation component in the den- 
sity distribution after last scattering. 



Paper III used the above methods to generate a model 
of a galaxy with a central black hole evolving from a small 
fluctuation at recombination. The final density profile 
was made a good fit to present day observations with data 
from M87, whose central black hole could be as massive 
as 3 X 10^ Mq. The initial fiuctuations at recombination 
were well within observational limitS"'^. Two possibilities 
were considered for the central black hole. The first sup- 
posed that it formed by collapse of matter during the 
evolution of the model, and for the particular model cho- 
sen, we found the black hole formed about 4 x 10® years 
before the present. The second supposed that it was a 
full Schwarzschild-Kruskal-Szekeres type wormhole with 
past and future singularities, but filled with matter [l^ . 
Again for the particular model chosen, it was found that 
the wormhole is only open for 6 x 10^^ s, and the original 
black hole mass was only 2 Mq, with all the rest of the 
mass accreting onto it over time. By recombination the 
mass was 2 x 10^ Mq, with a horizon 5 x 10^^ AU across, 
which corresponds to 4 x 10~^^ degrees on the CMB sky 
— much too small to leave an observable imprint (see 
Paper II for formulae relating the sizes of various objects 
to the angles that their images fill in the CMB sky). 



The subsequent sections complement these results by 
adding other criteria by which models can be specified, 
and deriving the solution algorithms that give the LT 
arbitrary functions in each case. 



^ The amplitudes of the density and of the velocity perturbation 
were within the limits set by observations done at the scale of 
f5! 1°. However, the scale appropriate for a single galaxy is 
0.004°, and at this scale there are no observational data at all. 



III. MODELS WITH A SIMULTANEOUS BANG 
TIME 

We here show how to find a LT model that evolves 
to a given 'final' (or 'initial' or 'middle') density profile 
starting from a simultaneous bang time 

t'g =Q, [ts = constant). (3.1) 

By (|2.2I) . the value of i? is unimportant near i? = 0, so the 
bang is Robertson- Walker (RW) like. This condition is 
known to generate only growing modes, but as discussed 
below, does not quite cover all possibilities, though the 
omitted case is not relevant to cosmology in an expanding 
universe. 



A. Density Profile Given at Time ti 

We specify a density profile Pi{M) at time ti, and 
choose tB- The function ai{M) is then determined from 
Pi{M) via (|2.24(l and H2.12(l . and the equations to be 
solved, in the three cases, are (ITTTj) . (|TTH|) & with 
ts — constant: 



HX: = \J [1 + aixY — 1 — arcosh(l + aix) 



dof 



x^/'^{U - tB) iPbdh{x); 



(3.2) 



EX: = arccos(l — aix) — ^1 — (1 — aixY 

-x^l'^it^-tB)"^ ^BDxix); (3.3) 



EC: — 27r — arccos(l — aix) + ^1 — (1 — a^a;)^ 



x^/^{t,~tB) =^ i^BDcix) 



(3.4) 



Here the age of the model, {ti — tB) is a free parameter 
that must be specified along with a^. Equations (|3.3|l - 
(|3.2|l may be solved numerically for x{M) and hence 
E{M). In addition, we give series expansions for so- 
lutions that are near to an expanding parabolic model 
(nPX), and near to maximum expansion in an elliptic 
model (nEM) at ti, because the above solutions would 
encounter numerical difficulties close to the parabolic and 
maximum expansion borderlines: 



nPX: 



nEM: 



20 1 


( drp 




3ai ' 


V TP 


+ 14) 


2 




indrf^ 


ai 


4 ' 

at 


25/2aiV2' 



where 



Ti = ti — tB, 
\/2 3/2 

rP = — a. , 
drp = Ti- Tp, 



Tn = TT 
dTN = Ti — TN- 



(I) 



3/2 



(3.5) 
(3.6) 

(3.7) 

(3.8) 
(3.9) 

(3.10) 
(3.11) 
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One may think of rp as the time it would take a paraboHc 
model to reach a^, and Tjq as the time to maximum ex- 
pansion if tti were the maximum a value. But note that 
dr^v is not the time since maximum expansion, because 
when drjv ^ 0, the model is not exactly at maximum 
expansion, so and rjv are less than their maximum 
values. Note too that H2.14|l is used to define x for the 
near parabolic series, so a negative x indicates a slightly 
elliptic model. 



1. Existence of Solutions 

We next consider existence for each M value point by 
point, as the solution type does not have to be the same 
at each point. Since the argument is very similar to that 
of previous papers, we will only here consider one case, 
and otherwise just summarise the conditions. 

We take {ti — ts) > 0, and, for the EX case we 
calculate 



V2at 



3/2 



HX: fe-ta) < 3' ; (3.15) 
PX: (i,-ts) = ^_^; (3.16) 
EX: ^^l^<it,^ts)<n[^f'; (3.17) 

(Q, \ 3/2 
2. 



EC: 



/a\3/2 



(3.18) 
(3.19) 



The parabolic and maximum expansion borderlines are 
sufficiently obvious that they will not be listed in subse- 
quent sets of existence conditions. 



B. Velocity Profile Given at Time ti 

For this scenario, the equations to be solved are H2.21|l . 
& with ts = constant: 



V'BDX(O) = 0, 



dip 



BDX 



dx 



BDX,x 



lim 4'BDX =7^-i2/a,f/^{t,-tB), (3.12) 

X — »2 / ai 

3/2 



— aix 2 



V'i3BX,x(0) = 0, 



lim tpBDX,: 



00. 



Clearly 'ipBDX,x is the product of a positive term {^/x ) 
and a monotonically increasing term (in square brackets), 
and therefore it can only change sign if the second term 
is negative at a; = 0, i.e. 



V2 



3/2 



< iU-tB). 



(3.13) 



Only in this case does tpBDX descend below before ris- 
ing to the limit H3.12|l . and only if this limit is positive, 
i.e. if 



HX: 



EX: 




— 1 — arcosh 



dof 



X-"\U-tB) = i^BVHix) 



bf +X 



b^ + x 



dcf 



X^/^iU - tB) =' i^BVxix) 



EC: = 27r — arccos 



x^'^iU - tB) i^Bvcix); 



dcf 



(3.20) 



(3.21) 




(3.22) 



where bi — {R^t)i/M^/^ , which are to be solved numeri- 
cally, and, from series expansions near the borderlines: 



nPX: 
nEM: 
where 



5&f fdrp 25dTj 



28r|, 



16 h, 64 b'i ' 
"""^'^+3^1+3^^'' 



(3.23) 
(3.24) 



3/2 



(3.14) 



does it create a root at x > 0. Otherwise ipBDX ei- 
ther rises monotonically from the intial at a; = 0, or 
it descends from and then rises to a negative limiting 
value. Therefore (|3.13|) and (|3.14|) are the necessary and 
sufficient conditions for a growing mode, still-expanding, 
elliptic solution to exist. The complete set of conditions 



Ti — ti — tB, 

4 

drp ^ n- TP, 

> 2/3 

I ^ 

Xjv = — 

■J- 2^/xW 
= — 



(3.25) 
(3.26) 
(3.27) 

(3.28) 
(3.29) 
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Here & is a kind of average velocity; if the model were 
exactly at maximum expansion after Ti , then a would be 
GN = 2/xjv and 6 = 



Notice that the model is only fully defined once K is given 
— in other words, the limiting late time RW model must 
be fully specified. 



1. Existence of Solutions 



The existence conditions for each case are: 

4 

36? 



HX: b,>0, {U-tB)>^\ 



EX: 6, >0, {U-tB)<^; 
EC: bi < 0. 



(3.30) 

(3.31) 
(3.32) 



IV. MODELS THAT BECOME 
HOMOGENEOUS AT LATE TIMES 

We next find a LT model that evolves from a given 
'initial' profile and approaches a RW model at late times, 
t oo. Only expanding hyperbolic (and parabolic) mod- 
els have an infinite future, and for these we require 

2E^K Af , K a constant. (4.1) 

Clearly, this case has only decaying modes. The time 
reverse of this case — the collapsing hyperbolic model — 
also satisfies the same condition but was RW like in the 
infinite past. See 's'l for examples. 

This condition may also be applied to elliptic models, 
and although the inhomogeneities due to tB.r 7^ do de- 
cay initially, other modes can grow as the big crunch is 
approached because tc.r 7^ 0. However, (|4.1|l has the ef- 
fect of making the lifetime T{r) along the dust worldlines 
a constant (see (|2.8|) '). Now, it is known that the set of 
necessary and sufficient conditions for no shell crossings 
in an elliptic model is {is,r/A^,r < 0,tc,r/Mr > 0} 
but T(r) = const makes it impossible to obey both these 
inequalities, unless Ib — const and tc — const, in which 
case the Friedmann model results. Thus the only elliptic 
model without shell crossings that satisfies 1)4. l|l is the 
RW modcl^. 



A. Density Profile Given at Time ti 
Applying l|4.1|l to 12.16|) leads to the direct solution 



HX: tB^U - y{l + a,KY-l 

-arcosh(l + a,if)]/i^3/2; 

nPX: tB^ti- ^[2^ 



a,K 3ajK^ 



20 



224 



(4.2) 
(4.3) 



^ Nevertheless, part of the evolution of non-homogeneous elliptic 
models obeying 14.11 will be free of shell crossings and may be 
of interest. 



1. Existence of Solutions 
The condition is 

HX: if > 0, (4.4) 

and it follows that 



\/2 3/2 

ti ~tB < ■ 



(4.5) 



This will keep appearing in what follows as a complemen- 
tary inequality to the others we shall derive. 



B. Velocity Profile Given at Time ti 

Given the velocity distribution hi{M), the hyperbolic 
case again has a direct solution for t^, once ()4.1|l is ap- 
plied. 



HX: 



is = ti 




arcosh 



nPX: 



Ib ~ ti 



4/1 2K 3K^ 



3 5bf 7bf 



^3/2; (4.6) 
(4.7) 



1. Existence of Solutions 



HX: >0, 0<K<bl {U-tB)>^. (4 



V. MODELS WITH A SIMULTANEOUS 
CRUNCH TIME 

For elliptic models, the requirement of only decaying 
modes is that the crunch time must be simultaneous (i.e. 
the crunch be RW like), 



tc — constant — tB 



2itM 
(-2£;)3/2 



= tB 



27r 



(5.1) 



H & P models with a simultaneous crunch time obviously 
have no bang, and are therefore collapsing at all times. 
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A. Density Profile Given at Time ti 

The equations to be solved are l|2.17|) & (|2.18|l with 
l|5.HI and tc = constant: 

EX: — —2tt + arccos(l — aix) — ^^1 — (1 — aixY 



dcf 



EC: 



arccos(l — a^x) + -^/l — (1 — a^x)^ 



dcf 



HC: = v^l + OiXp — 1 — arcosh(l + aix) 



dcf 



X^/'^{tc - U) ljjCDH{x) 



with the equations for the borderhne cases being: 



(5.2) 



(5.3) 



(5.4) 



nPC: xK — i - 



nEM: x : 



20 
2 



25dT|, 



where 



n = tc — ti, 



(5.5) 
(5.6) 



(5.7) 



while (|3.8|I - H3.11|I define rp, tn and drAr. Again the re- 
maining lifetime of the model, {tc — ti) is a free param- 
eter that must be specified along with a^, and as before, 
eqs. (|5.2|l - (|5.4|l are to be solved numerically for x. The 
functions E and ts then follow from 12.14|l or H2.13|l and 



Existence of Solutions 



The complete set of conditions is 



EX: 7r(^j <{tc~t,) 

{tc 
V2 



EC: ^a,'/-'<{tc-t,)<n[^]'^' 



HC: {tc - *0 < — a. 



3/2 



(5.8) 
(5.9) 
(5.10) 



The first inequality in H5.9|l means that the time- 
difference between ti and tc is larger than it would be in 
a collapsing parabolic model — this is seen from H2.5|l . 
If {tc — ti) is smaller, then the model collapsing from 
the given density distribution to the Big Crunch in the 
time-interval {tc — ti) would have to be hyperbolic. This 
is consistent with the time-reverse of Eq. (|4.5|) . 



B. Velocity Profile Given at Time ti 

These models are found by solving 



EX: = —2tt + arccos 



+ x^^^{tc~t,) iJcvx{x) 



+x 



dcf 



-\ 1 



bf+X, 

(5.11) 



EC: 



HC: = 



6? 



bl 



6? 



dcf 



-| - 1 - arcosh -5 

bj-xj \bf-x 

dcf 



-x''/^{tc-t,) =i^cvH{x); 

or, for the borderline cases, by calculating 

5^2 fdrp 25dr|,^ 
nPX: X ' 



6 \tp 28t|, J ' 

,1 16 64 6f 

nEM: - - ( 1 + + ^=1 ; , 



where 



Ti = tc — ti, 



(5.12) 
(5.13) 

(5.14) 
(5.15) 

(5.16) 



while IpO^ll and define rp, xn and 6. 



i. Existence of Solutions 



EX: 6, > 0, tc> tf, 
EC: 5, < 0, {tc - t,) < 



HC: b, < 0, (tc - t,) > 



3br 
4 



(5.17) 
(5.18) 

(5.19) 



where the last equation above is for collapsing hyperbolic 
models. By ^^,x^{tc - ti) = 2-n - x^l'^{ti - ts), 
and using this in (|5.11|) leads to an equivalent existence 
condition for the EX case, 6; > 0, (ti — tg) > 

This means the conditions (|4.8|l and H5.17|l are mutually 
exclusive. 



VI. GROWING AND DECAYING MODES 

For pure decaying modes, hyperbolic (and parabolic) 
models that are expanding must become RW like at late 
times, whereas for elliptic models and for hyperbolic and 
parabolic models that are collapsing, they must become 
RW like at the crunch. 
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Conversely, for pure growing modes, hyperbolic models 
that are collapsing must have been RW like in the distant 
past, while elliptic models and expanding hyperbolic and 
parabolic models must have been RW like at the bang. 

This is summarised in the following table: 





only 


only 




growing 


decaying 




modes 


modes 


HX, PX 


tB,r ^ 




E 


tB,r = 


tc,r = 


PC, HC 


2E = KM"^!^ 


tc,r = 



VII. 



MODELS WITH A SIMULTANEOUS TIME 
OF MAXIMUM EXPANSION 



Not infrequently, authors seeking a manifestly regular 
initial condition for an inhomogeneous matter distribu- 
tion, have required a finite density distribution and zero 
initial velocity. This is achieved in an LT model if the 
moment of maximum expansion occurs at the same time 
along all worldlines. 

Naturally only elliptic (EX & EC) models have a mo- 
ment of maximum expansion. In general, the time of 
maximum expansion 12.9|l is different for each worldline. 
The condition that it be simultaneous is ijv/x = constant, 
i.e. 



tsMX = constant = 



(-2£;)3/^ 



= tB 



j.3/2- 



(7.1) 



This needs to be combined with another condition to 
obtain a solution. 



B. Density Profile Given at a Simultaneous Time 
of Maximum Exapnsion 



Although this is a special case of section IVII Al it has 
an explicit solution. Given p = psMx{M) ai t ~ tsMX, 
we calculate Rsmx{M) with H2.24|l . then write (|2.1U|) as 



. , M 2 
{-E) = ^ X : 



RsMX asMX 

and use it in IjT.lfl . giving the direct solution 



is = tsMX — TI" 
M2/3 

E — . 



'SMX 



8 



asMX 

Solutions obviously exist for Rsmx > and M > 0. 

C. Velocity Profile Given at Time ti 

The two elliptic cases are found from. 



(7.6) 

(7.7) 
(7.8) 



EX: = — TT + arccos 



+ x^/^{tsMx - U) iJsvx{x)- (7.9) 



5? 



b'? + x 



dcf 



EC: = TT — arccos 



+ x^/\tsMX - U) = i^svcix). (7.10) 




1 - 



bj-x_ 
bj +x 



dot 



A. Density Profile Given at Time ti 

The equations to be solved, with tsMX — constant, are 



EX: = — TT + arccos(l — aix) — ^^1 — (1 — aixY 

+ x'^'^itsMX - U) =' i^sDx{x)- (7.2) 
EC: = TT — arccos(l — a^x) + ^1 — (1 — aix)'^ 

+ X^'/^itsMX - U) =' ^SDC{X). (7.3) 

As above, these equations are solved numerically for x, 
and E and t b then follow from H2.13|l and (|7.1|) . 



1. Existence of Solutions 

The complete set of conditions is 

EX: b, > 0, tsMX > U; (7.11) 
EC: b, < 0, tsMX < U. (7.12) 

The conditions 6^ > and 5^ < follow directly from the 
assumptions that the model is expanding or collapsing, 
respectively. The other ones are conditions for solvability 
of the corresponding equations. 

VIII. MODELS WITH GIVEN DENSITY AND 
VELOCITY PROFILES AT THE SAME TIME 



1. Existence of Solutions 



The complete set of conditions is 



EX: tsMX > U; 
EC: tsMX < U. 



(7.4) 
(7.5) 



For this case, we are given R^t{M,ti) and p{M,ti) 
which provide us with bi{M) and ai{M) via equations 
and (|TT5|l . We then solve for E, 



2E^iR,t)^-^ ^ ±x^ 
Hi 



2E 

M2/3 



b 



(8.1) 
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and obtain from one of (|2.16() - (|2.19|) . The equation 
for tB is sensitive to the sign of bi (i.e, i?,t), even though 
the E equation isn't. This method is given in [T?il |. 

The simple solution H8.1() does also follow as the limit 
t2 — > ti from the considerations of Paper II, though 
merely substituting t2 = ii leads to degenerate equa- 
tions. See Appendix A for a proof. 

The borderline cases need no special treatment, as 
there are no numerical difhculties arising from being close 
to them. The model is parabolic if 6,f — 2/ai, and is at 
maximum expansion (at ti) if bi = 0. 



A. Density Profile Given at Time ti 

Since x is known from H9.,S|I and ti, pi and Oi are finte, 
we find tB from: 



HX: tB^t,-[^{l + arblJ^-l 

- arcosh(l + a,6Le)] /^fate; (9.4) 
nPX: tB^t,-Tp{l-^^^^+^-^^^- (9.5) 



20 224 J 

which are well defined for all 6iate > 0, but imply (|4.5|l 



1. Existence of Solutions 

Eq H8.1|l always has a solution, but, for the purposes 
of determining ts, the various types of solution exist if 



HX: 


h^> 


2 


and \ 


li > 0; 


(8.2) 






fli 




EX: 


bi< 


2 


and \ 


h > 0; 


(8.3) 






ai 




EC: 


h\< 


2 


and \ 


li < 0; 


(8.4) 






ai 




HC: 




2 


and \ 


h < 0. 


(8.5) 






ai 







1. Existence of Solutions 
The model is expanding if 

hate > 0, 

and the above solution is well defined provided 

> 0. 

B. Velocity Profile Given at Time ti 

The equation for is 



(9.6) 



■7) 



IX. MODELS WITH A GIVEN VELOCITY 
PROFILE AT LATE TIMES 



By "late times" we mean the asymptotic future, i.e. 
the limit ry — > oo and t ^ oo, so this section applies 
only to expanding hyperbolic (HX) models. The time 
reverse — a collapsing hyperbolic (HC) model with a 
given density or velocity profile in the infinite past — 
follows in an obvious way. 

From lEini) and we find 



R 



lim 

77^00 t — tB 

lim R t — 

t — >oo ' 



(9.1) 
(9.2) 



which gives simply 



HX: tB — ti — 

— arcosh 
nPX: tB^U-Tp(l 



^fatc 



- 1 



56? 



"late 
_ ^"latc 

7bf 



"late I 



(9.9) 



1 . Existence of Solutions 

This solution is well defined (for an expanding model) 
provided 



bi > 6iate > 



(9.10) 



This does mean that not all choices of 6iatc(Af) and bi{M) 
are possible, but it should be easy to choose a suitable hi 
after 6iate is fixed. 



E 



R.t 



^ ^ ^ - &latc- 



(9.3) 



C. Simultaneous Bang Time Specified 



This always exists, and fully determines E{M), leaving 
tB{M) free to be determined by a second requirement. 



As soon as the (constant) value of is given, we have 
all the LT arbitrary functions E{M) and tB{M). 
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1. Existence of Solutions 
Existence only requires that the model be expanding, 

6late > 0. (9.11) 



X. MODELS WITH A GIVEN DENSITY 
PROFILE AT LATE TIMES 

From (EH) and wc obtain, using dR/dM = 



dR/dM 
{t - ts) 



'2E 



dE/dM\ (cosh 77- 1) 
E J (sinh rj — rf) 



1 

M 

3dE/dM 1 \ sinh?7 
2E MJ (cosh 77- 1) 

{2EfdtB/dM sinh 77 



M 
2 



/ dR/dM \ 



(coshTy — 1) (sinh 77 — 77) 
(10.1) 

(10.2) 



^t-tB) J \ (t-ts) J 

Assuming dt b /dM is finite^ , we find that 

dR/dM dE/dM 
lim = — , , 

'V^oo [t-tB) y/2E 

which with 19.1|l leads to 
lim [Kp{t - ts)^] = 



y/^i:E dE/dM 



(10.3) 



(10.4) 



Clearly p(t — is)^ freezes out — becomes time indepen- 
dent — and 



lim Kp{t — tB}'' 

rj — >oo 



6 



((2£;)3/2))' 



(10.5) 



dAI 



Thus, if we specify the late time limit [p(t — tB)^]\atc{M), 
we have 



3/2 _ {2E) 



3/2 



''late 



M 



M 



6 



1 



M Jq K[p{t ~ tB)^]utc 



dM. (10.6) 



This assumes an origin exists, i.e. £^ = at Af = 0. If 
not, then we must specify some E — Ei aX some M = Mi. 

Again, (fTITHjl fully determines E{M), leaving tB{M) 
free. 



^ One can consider models in which dtg/dM becomes infinite while 
dE/dM doesn't, either asymptotically or at individual points. 
However, at such locations there is a permanent zero in the den- 
sity, and the above limit is not valid. The other cases that give 
a permanent zero in the density, loci where dE/dM is divergent 
and dtg/dM isn't, and loci where M,r = are not problematic 
for this limit. (See Paper II for a discussion of regions of zero 
density. ) 



A. Density Profile Given at Time ti 

In this case we find Ib from: 



HX: tB - U - [V(l + a»a^iato)2-l 

- arcosh(l + aixutc)] / xfl^^; (10-7) 



nPX: tB^U-Tp{l- 



20 



224 



;(io.8) 



which is well defined for all xiatc > 0, but again implies 
(lOll . 



1. Existence of Solutions 

Equation (|10.6|) will have a positive solution for x for 
any positive [p{t — iB)'^]iatc(-^) and the above equation 
for is well defined for any aj that is derived from a 
positive Pi. 



B. Velocity Profile Given at Time ti 

The equation for is 



HX: tB — ti 



-^latc 



1 



bi — ablate , 

bf + Xlato \ 



arcosh 



•^late y 



nPX: t 



B 



ti-Tp [1 



9x 



late 



5&? 



7bj 



3/2 , 

''late' 

(10.9) 
(10.10) 



1. Existence of Solutions 

The above is well defined provided 
> a:^iato- 



(10.11) 



It should not be hard to choose bi to satisfy this, once 
xiatc is known. 



C. Simultaneous Bang Time Specified 

Again, setting the constant value of completes the 
specification of the LT arbitrary functions, and the solu- 
tion always exists. 



D. Late-Time Velocity Profile Given 

Since specifying the late-time velocity profile via eq. 
(|9.3|) . and specifying the late-time density profile via 



11 



eq. 
ble. 



(|10.6() both fix a;, this combination is not possi- 
(Should it happen that the two specifications are 
consistent, then another 'boundary condition' would be 
required to specify a particular LT model.) 



XI. CONCLUSIONS 



We have developed several new ways of specifying the 
'boundary' data needed to uniquely determine the evo- 
lution of a Lemaitre-Tolman model, and thus provided 
more options for designing models with particular prop- 
erties or behaviours. Thus one can now easily generate 
models that start from an initial stationary state, or have 
only growing modes, or approach a specified density or 
velocity profile in the asymptotic future, or approach RW 
models at late times or diverge from them, etc, as listed in 
the introduction. The foregoing properties are combined 
in pairs to fully specify a particular LT model. Although 
several of the individual properties considered here have 
previously been used, what is significant here is the equa- 
tions that result from the many combinations of pairs of 
properties, and the derivation of existence conditions for 
the 3 types of solution. Also, our results have been pre- 
sented in a form that is easily converted into coding for 
numerical calculations. Current work [T^ provides an 
example of the use of some of the new LT specification 
methods to create and evolve a model of the Shapley 
concentration and the great attractor. 
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APPENDIX A: DERIVATION OF FROM 
PAPER II RESULTS 



Using the results of Sec. IV of Paper II with t2 = ti = 
ti, and remembering that for the EC case both profiles 



must be specified in the collapsing phase, we obtain 



HX: 




= Y (1 + aix) — 1 — arcosh(l -|- aix) 

'bf +x 



cosh 



bf-x^ 
EX and EC: 

= arccos(l — aix) — ■\/l — (1 — Uix)^ 



(Al) 



+ yi- 

Denote 1 -f- aiX 

hyperbolic case, 
sinhit— u — sinh 
is single-valued, 
one solution, u 
bi^ — 2/ai + X. 
case, this result 
The result bi" 
X = -2£;/M2/3, 



bf +x 



— arccos 



b'f + x 



(A2) 



= coshu and /'i — coshw for the 
Equation (|AHl then is equivalent to 

v—v. But the function := sinhy—y 

so the equation F{u) — F{v) has only 
= w. In our case, this solution is 
bmce X — 2E/M'^/'^ in the hyperbolic 

is equivalent to H8.1|l . 

= 2/ai — X for the elliptic cases, where 

follows in a similar way. 



APPENDIX B: CALCULATION OF SERIES 
EXPANSIONS 

Most of the series expansions given above are non- 
trivial to derive, so a couple of representative calculations 
are outlined here. 

In each case, such as equations (|3.2f) - (l3.4|) . it is usu- 
ally possible to do a direct series expansion in powers of 
some suitable small variable, and then invert it to get a 
series for x. However, the result is usually not very tidy, 
and tends to have sign ambiguities. Much better defined 
and neater results always follow from starting with the 
parametric evolution equations, H2.4|l - H2.6|l . The calcula- 
tions were all done to 6th order, using Maple, but only 
truncated series are written out here. 

Let us obtain the near-Parabolic series that lies be- 
tween (|3.2ll and 1)3. 3|l . The parabolic limit occurs when 
a; — > 0, while R and t = t — ts remain finite. This 
requires 



rj ^ and 



V 



(Bl) 



so that the new evolution parameter e remains finite for 
finite T. Series expansions of (|2.()|l then give 



(B2) 
(B3) 



e3 


xe^ 
1 


x^e? 


X e: 


6 ' 


120 


^ 5040 


^ 362880 




4 

xe 


' 1 


x^e^ 


2 




720 


40320 ^ 



Now we invert the above series for r by writing 



eo -I- e\x -I- 623; 



(B4) 
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substituting into ljB2|) . and solving each power of x in 
turn for the cocfhcients e^. We get 



(6t)4/3x2 (6T)2a;3 



1400 25200 
which we substitute into (jBSp . 

2 20 

3(6r)4/3a;2 23(6r)2a;3 



(B5) 



2800 



504000 



(B6) 



We define rp as the time (since the bang) that it would 
take an exactly parabolic model to expand to a 



{6tp) 



(B7) 



and we define drp as the difference 

dr = T- Tp. (B8) 
Then we once more invert the series by writing 

X « xidr + X2dT'^ + x^dr^ + ■ ■ • , (B9) 



substituting into (|B6p . and solving each power of dr in 
turn for the coefficients Xi, thus obtaining 



20 / drp 2bdTl lOOOOdr^, 



3a 



Tp Urf 



3969r3 



, (BIO) 



which is H3.5|l . 

As our second example we find the near maximum ex- 
pansion series lying between lpOT|l and (j2!22l- Maxi- 
mum expansion occurs when 77 — > tt, 6 — > 0, a — > 2/x, 



and r tt/x^/^. The mean velocity is 



2 / TT 



1/3 



Thus we write 

T] ^ n + e, 
expand in powers of this new e, 



4 24 



and invert to get 



\px \ 3x 



(Bll) 

(B12) 

(B13) 
(B14) 



(B15) 



We then substitute this as well as 

xK,xxh + xih^ -I- xsfo^ H , (B16) 

into (|B13|) . to get the second series inversion: 

20 / 166 6462 I28(97r2 - 112)6^ \ 

XK. — 1 H = -I ^ i = '- \ , 

(B17) 

where 



2 /vr 



TT \ Tj 



1/3 



(BIS) 



is what the average velocity would be if maximum ex- 
pansion occurred at r^. This gives us H3.24|l . 
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